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ABSTRACT 

We investigate the idea that the interaction of Dark Matter (DM) sub-halos with the gaseous disk 
of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) 
distributions. We use high resolution hydrodynamic simulations to show that pure DM sub-halos 
impacting a galactic disk are not able to produce holes; on the contrary, they result in high density 
regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total 
DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The 
size and lifetime of these holes depend on the sub-halo gas mass, density and impact velocity. A DM 
sub-halo of mass 10* Mq and a gas mass fraction of ~ 3%, is able to create a kpc scale hole, with a 
lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation 
rate at the rim of the hole, again in agreement with observations. Even though the properties off these 
simulated structures resemble those found in observations we find that the number of predicted holes 
(based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) 
falls short compared to the observations. Only a handful of holes are produced per Giga year. This 
leads us to conclude that DM halo impact is not the major channel through which these holes are 
formed. 
Subject headings: galaxies: structure — ISM:HI — ISM:bubbles — cosmology: dark matter 



1. INTRODUCTION 
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distribution of many nearby 
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__ ._ ... MlOl and NGC 6946 

(|Kamphuislll993D. IC lOdWilcots fc Milled [l998h , SMC 
(|Stavelev-Smith et al.l 11997^ ■ LMC (|Kim et al.l [T9M) 
and IC 2574 ([Walter fc Brinks! [19991 ) . These structures 



are typically regions of low HI density distributed across 
the entire galaxy, and come in various sizes ranging from 
a few 100 pc to about 1.5 kpc. Their expansion velocity 
sets an upp er limit to their dynami cal age of about 10 - 
60 Myr (cf. iWalter fc BrmS (fT999h for IC 2574). 

The formation of HI holes and shells has often 
been ascribed to the action of stellar winds and 
supernova explosions occurring in OB asso c iation s 
and young stellar clusters (c f. iWeaver et al.l (I1977D. 
iTenorio-Tagle fc Bodenheimeii ()1988[) . Ivan der Hulstl 
(J1996D ). Bas ed on a simple mode l of holes created by O 
and B stars, lOev fc Clarkd (J1997D successfully predicted 
the observed number distribution of holes in the SMC, 

lending support to thi s hypothes is. Mor e recently, the 

analys is conducted by lOtt et al.l (I2001D. iSimpson et al.l 
(l2005l) . iWeisz et al.l (|2009aD and ICannon et al.l (|2011D 
have shown that the stellar content of HI holes in 
Holmberg I and II, DDO 88 and 165 can release enough 
mechanical energy to drive the formation of holes on 
timescales similar to what is observed. In particular, 
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iWeisz et al.l (|2009a[ ) were able to show that all the holes 
in Holmberg II observed with H ST contain multiple 
stellar populations of different age. I Warren et al.l (|2011[ ) 
studied the HI holes in five dwarf galaxies (DDO 181, 
Holmberg I, M81 Dwarf A, Sextans A and UGC 8508), 
and failed to detect a young star cluster at the center 
of the observed HI holes. They thus suggested that 
large holes may form due to multiple episodes of star 
formation. 

There is more observational evidence against a single 
burst of star formation being responsible for the creation 
of HI holes. F or example , no ro bust spatial correlation 
was found by iKim et aLl (|1999f ) between the distribu- 
tion of the Ha emiss i on an d the HI holes in the LMC. 
iHatzidimitriou et al.l (|2005[ ) estimated that the holes in 
the SMC not associated with a young star cluster are a 
factor of 1.5 more than those exhibiting relatively young 
stars at their center. In the Galaxy, the distribution of 
radio pulsars compared to that of holes seem to indicate 
that holes can not be produced by th e supernova explo- 
sions o f a single age stellar po pulation (jPerna fc Gaenslerl 
I2004D . iPasauah et~all (I2008D found that the most recent 
episode of star formation has taken place preferentially 
at the rims of the HI holes in IC 2574, exception made 
for the supergiant shell which indeed emb eds a young 
star cluster fcf. iStewart fc Waited (|2000[) . IWeisz et al] 



(l2009b[) V iRhode et all (|1999[) observed the HI holes in 
Holmberg II in search of the descendants (stars of spec- 
tral types B, A and F) of those clusters whose supernovae 
would have produced the holes. They measured an inte- 
grated light of the stars within the holes inconsistent with 
the hypothesis of a young star cluster triggering the for- 
mation of a HI hole. In the case of IC 1613, its largest HI 
shell was found to host about 27 OB associations whose 
energy ca n not sustain the formation and expansio n of 
the shell (jBorissova et al.l (|2004[) . ISilich et al.l (|2006[ )). 
Some authors have proposed alternative formation 
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hyp otheses to the SNe origin. lEfremov et al.l (|1998[ ) 
and iLoeb fc Pernal ()1998D postulate that a high-energy 
gamma ray burst (GRB) from the death of a single 
massive star could create kpc sized holes in the Inter 
Stellar Medium (ISM), thus offering an explanation for 
holes without a detectable underlying cluster. These au- 
thors assume the energy from GRBs is emitted isotropi- 
cally. However, GRBs r elease most of their energy in bi- 
directional beams (eg., iBlandford fc Znaie3ll977l mech- 
anism), making this scenario less likely to produce large 
HI holes. 

Another mechanism prop osed is the infall of gas clouds 
(|Tenorio-Tagle et al.l 119871 ). One observational predic- 
tion of this model is a half-circle arc seen in an HI 
position- velocity diagram. The half-circle arc arises from 
the gas being pushed to one direction, corresponding 
to the direction of travel of the high velocity clouds. 
Some o bservational s upport for th is idea is reported by 
Heiles (jHeilesI 119791 iHeilesI I1984D who point out that 
the most energetic Galactic shells in their study all 
h ave half-circle arc sign atures in position-velocity space 
( IKamphuiset al.lll99ll ). 

In this paper we investigate an alternative formation 
scenario, where holes and shells in extended HI disks arc 
the result of interaction of the gaseous disk with dark 
matter (DM) sub-halos. The Lambda Cold Dark Mat- 
ter (LCDM) model predicts the existence of thousands of 
dark matter substruc tures within the dar k matter halo 
of every galaxy (e.g. ISpringel et al.l [20081 and reference 
therein). The majority of these sub-halos is "dark" i.e. 
does not host a satellite galaxy or any visible ste llar 
structures (e.g. iMaccio et al]|2010l Hb^t et al.ll20lH) . If 
these sub-halo population is able to produce the observed 
holes in galactic HI disks, like the one of IC 2574, then 
this might provide a new way of detecting the presence 
of non-luminous DM halos. For example, the number 
of holes could then be linked statistically to the amount 
of substructure in the dark matter distribution and thus 
provide a unique way to asses the nature of dark matter 
and by extension to test the ACDM model. 

Perturbations on a stellar disk due to (massive) satel- 
lites have been extensively studied in the recent y ears 
(e.g. iKazantzidis et al.l 120081 iMoster et all l2010al and 
references therein) , while less attention has been given to 
the effect of low mass DM clumps on an extended gaseous 
disk, with few exceptions. iBekki fc Chibal (|2006( ) investi- 
gate how the impact of dark matter sub-halos orbiting a 
gas-rich disk galaxy embedded in a massive dark matter 
halo influences the dynamical evolution of the outer HI 
gas disk of the galaxy. They show that the impact of dark 
matter sub-halos ("dark impact") can be important for 
better understanding the origin of star formation discov- 
ered in the extreme outer regions of disk galaxies. They 
also discuss the possibility that this kind of dark impact 
will be able to produce holes in the gas distribution. In 
their study they adopted a model, which did not include 
multiple events, a live dark matter halo (they adopted 
an analytic fixed potential), cooling and star formation. 

A new attempt to detect the imprint of the dark satel- 
lites in the HI disk of the Milky Way has been r e cently 
mad e. In a series of papers iChakrabarti fc Blitd ()2009[ ) 
and IChakrabarti fc Blitd ()2011t) examine tidal interac- 
tions between perturbing dark sub-halos and the gas disk 
of the Milky Way using high-resolution Smoothed Par- 



ticle Hydrodynamics simulations. They compare their 
numerical results to the observed HI map of the Milky 
Way, and find that the Fourier amplitudes of the pla- 
nar disturbances are best fitted by a perturbing dark 
sub-halo with a mass that is one-hundredth of the Milky 
Way with a pericentric distan ce of 5 kpc. More recently 
( iChang fc Chakrabartill20Tll ) develop a perturbative ap- 
proach to study the excitation of disturbance in the ex- 
tended atomic hydrogen discs of galaxies produced by 
passing dark matter sub-halos. They show that the prop- 
erties of dark matter sub-halos can be inferred from the 
profile and amplitude of the different perturbed energy 
modes of the disk. 

Motivated by these recent studies, in this paper we use 
high resolution hydro-dynamical simulations to study in 
details the interaction of DM sub-halos with a galactic 
gaseous disk. Our primary goal is to see under which 
conditions DM satellites are able to create holes that re- 
semble the observed ones and whether the majority of 
these holes can be explained by dark satellites-disk in- 
teractions. We model our primary galaxy on the nearby 
dwarf galaxy IC 2574 and try to replicate its observed 
features in our simulations. After having described our 
numerical setup, we will first present results of a single 
satellite-disk interaction, for different satellite orbital pa- 
rameters and gas content. Then we will use satellite pa- 
rameters (mass, velocity and position) directly obtained 
from high resolution N-body simulations to study the fre- 
quency of DM sub-halo disk encounters. Finally we will 
discuss the implications of our results. 

2. NUMERICAL SIMULATIONS 

We make use of the parallel TreeSPH-code GADGET- 

2 (Springer "2005^ in this work T he cod e uses SPH 

(iGingold & Monaghan 19771; ILucvl 119771: l^lonagha^ 



Il992f ) t o evolve the gas using an entropy conserving 
scheme (jSpringel fc Hernguistl 1200^ . Radiative cooling 
is implemented f or a primord i al mix ture of hydrogen and 
helium following iKatz et all (I1996D and a spatially uni- 
form time-independent l ocal ultraviolet backgro und in 
the optically thin limit (jHaardt fc Madaul 119960 is in- 
cluded. The SPH properties of the gas particles are aver- 
aged over the standard GADGET-2 kernel using 64 SPH 
particles. Additionally, the minimum SPH smoothing 
length is required to be equal to the gravitational soft- 
ening length in order to prevent artifi cial sta bilization 
of sm all gas clumps at low resolution (jBate fc Burker^ 
\T99Th . 

All simulations have been performed with a high 
force accuracy of aforce — 0.005 and a time intcgra- 
tion accuracy of rjacc = 0.02 (for further details sec 
ISpringelll2005[ ). Star Formation (SF) and the associated 
heating by supcrnovae (SN) are modeled following the 
sub-resolution mult iphase interstellar m e dium (ISM) 
model described in ISpringel fc Hernguistl ()2003[ ). The 
ISM in the model is treated as a two-phase medium 
with cold clouds embedded in a hot component at 
pressure equilibrium. Cold clouds form stars in dense 
(p > pth) r egions on a tim e-scale chosen to match 
observations (jKennicuttI I1998D . The threshold density 
Pth is determined self-consistently by demanding that 
the equation of state (EOS) is continuous at the onset 
of SF. Wc do not include feedback from accreting black 
holes (AGN feedback) in our simulations as there is 
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no evidence of AGN activity in dwarf galaxies like the 
one studied here. In our runs, the parameters for the 
star formation and feedback are adjusted to match 
the Kroupa initial mass function (IMF) as specified by 
iKroupal (|2001f) . The SF time-scale is set to i^ = 3.5 
Gyr, the cloud evaporation parameter to Aq = 1250 and 
the SN temperature to Tsn = 1-25 x lO^K. 



2.1. Primary Galaxy Setup 

We apply the method given by iSpringel et al.l (|2005l ) 
to construct the central galaxy. The central galaxy con- 
sists of an exponential stellar disk. In order to match 
the almost constant radial gas density profile of IC 2574, 
we have modeled the gas disk by two components. A 
radial exponential component and a constant radial pro- 
file which falls off rapidly at a specified scale radius. The 
stellar disk has a mass Mdisk, the gaseous disk has a mass 
Mgas, with a spherical bulge with mass Mb em bedded in a 
dark m atter halo of mass M^h-. The halo has a lHernguistl 
(|1990( ) profile with a scale radius corresponding to a 
Navarro Frenk & White halo (NFW: [Navarro et al.lll997| ) 
with a scale-length of rg and a concentration parameter 
Cvir = Rvir/rs- Wc usc the results of^laccio et al. (2008J) 
to compute halo concentration as a function of virial 
mass. The scale-lengths r^ of the exponential gaseous 
and stellar disks are assu med to be equa.1, and are deter- 
mined using the model of lMaccio et al.l ()2008[ ). assuming 
that the fractional angular momentum of the total disk 
Jd = (Jgas + ^disk)/>/vir is cqual to the global disk mass 
fraction m^ — (Mgas + Mdisk)/Mvir for a constant halo 
spin A. This is equivalent to assuming that the specific 
angular momentum of the material that forms the disk 
is the same as that of the initial dark matter halo, and 
is conserved during the process of disk formation. The 
vertical structure of the stellar disk is described by a ra- 
dially independent sech? profile with a scale-height zq j 
and the vertical velocity dispersion is set equal to the 
radial velocity dispersion. The vertical structure of the 
gaseous disk is computed self-consistently as a function of 
the surface density by requiring a balance of the galactic 
potential and the pressure give n by the EOS. Th e stellar 
bulge is constructed using the iHernguis^ ()199C1[) profile 
with a scale- length r^ . 

We build a primary galaxy which match es the prop- 
erties of IC 2574 (mass profiles taken from iLerov et al.l 
I2008J) . The galaxy has a halo of mass M200 = 10" M© 
obtained from the st ellar mass, from the recipe given in 
iMoster et al.l ()2010bl ). containing a stellar disk of mass 
Mdisk = 3.16 X 1O®M0, a gaseous disk component of mass 
Mgas = 1-8 X IO^Mq with 80% of the gas in the HI disk, 
a small bulge of mass Mbuige = 3.16 x IO^Mq. The DM 
halo has a concentration parameter of C200 = 8.47. The 
scale radius of the disk (stellar & the exponential com- 
ponent of the gaseous disk) is rd = 2.1 kpc and the scale 
radius of the slab of gas is Tsiab = 4.5 kpc. Once we fix 
rd we compute th e halo spin parameter from the recipe 
of iMo et all (|1998[ ) which resuhs in a value of A = 0.047. 
The disk scale- height is fixed at zq = 0.15 kpc and scale 
radius of the bulge is set at rb = 0.095 kpc. The galaxy 
has iVDM = 10^ DM particles, iVdisc = 10^ in the stel- 
lar disc, A'^buigo = 10* in the bulge, and A'^gas = 5.7 x 10^ 
gaseous disk particles. The force resolution (softening) is 
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Fig. 1. — The red line is the gas surface density and the black Hne 
is the surface density of stars in the inner 7 k pc of IC 2574. Th e 
observational data is taken from Figure 27 of l|Lerov et al.ll2008l V 
The error bars represent the rins uncertainty value. 



101,80,57 pc for dark, gas, and stars respectively In order 
to have a stable initial condition we evolve this primary 
galaxy in isolation for 1 Gyr. The initial conditions were 
chosen in such a way that the surface density of the gas 
match the observed density of IC 2574 after evolving for 
1 Gyr as shown in Figure [T] 

3. SIMULATION RESULTS 

We here present the results of our numerical experi- 
ments. We mainly run three kinds of simulations: single 
encounter with a pure DM sub- halo (no gas or stars); 
single encounter with a DM sub-halo containing a small 
gas fraction; and multiple, cosmologically motivated en- 
counters. 

These single halo simulations are simplistic experi- 
ments designed to explore the effect of sub-halo/disk in- 
teractions. Orbit, mass and gas content are not meant 
to represent the typical case. An cosmological motivated 
run with multiple interactions will be presented in sec- 
tion [3?3l 

3.1. Pure DM Sub-halo Interaction 

In this section we investigate the dynamical impact of 
a pure DM halo on the gaseous disk. The mass of the 
sub-halo is first fixed at a value of Afgb = IO^Mq. We 
want to have a highly concentrated sub-halo to act like 
a 'bullet' in the interaction. To do this we start with a 
^200 = 10^°M(3 halo and then carve out 99% of its mass 
and create a 10 M0 halo. This is justified, as it has been 
shown that the sub-halos passing through the primary 
halo of a galaxy get tidally stripped and can lose up to 
90% and in some e xtreme cases 99% of their mass dur- 
ing t heir orbit (see iPeiiarrubia et al.ll200"8t iMaccio et al.l 
l2010f ). The total number of DM particles in halo has been 
set to Ndm = 1000. Starting from a lower mass halo (e.g 
IO^Mq) would have resulted in a less concentrated 'bul- 
let', hence our choice maximizes the dynamical effect on 
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Fig. 2. — Surface density map of the gaseous disk after the pas- 
sage of the pure DM sub- halo of mass 10* M©. Such an encounter 
leads to a gas overdensity (marked by a circle) . The units of den- 
sity in this plot are in lO^'^ MQkpc~^ . The 2D contour plots were 
made using the visualization software SPLASH (Price 2003). 



the disk. This sub-halo is placed at a distance of 5 kpc 
above the gaseous disk of the primary galaxy and given 
an initial velocity of w^ = 150 kms~^ perpendicular and 
pointing towards the galactic disk. Figure [21 shows a sur- 
face density map of the gaseous disk after the passage of 
the sub-halo. 

As it emerges clearly from Figure [2J a pure DM sub- 
halo is unable to form a hole, instead it gives rise to a 
localized high density region (marked by a circle). This 
is due to the fact that the DM particles are collision less. 
Due to the lack of contact forces the sub-halo cannot 
push away the gas particles. DM can only gravitation- 
ally focus the gas in a stream behind its path of motion 
(Bondi-Hoy le accretioniEdgar 2004: Bondi &: HovldllQll 
IBondi|[r952l ) which gives rise to higher surface density in 
the region it passes through. These high density peaks 
are 20% — 25% denser than the mean density of the disk, 
and they last for about 80 — 90 Myr, after wards they are 
destroyed by the dispersion of the gas and the differen- 
tial rotation of the disc. The fundamental result of this 
first experiment is that a pure dark matter sub-halo is 
not able to create a hole in the gaseous disk. 



3.2. Go 



Sub-Halo Interaction 



Since a pure DM halo is not able to displace the gas 
in the disk, a medium which interacts through gravity 
as well as contact forces is needed in order to push the 
gas away from its motion path. A gaseous medium will 
provide the required contact forces to push the gas away 
and form holes. We construct the DM sub-halo with the 
same mass and concentration parameter as described in 
the previous section and in addition we also add a hot 
gaseous component. The hot gas has a beta profile with 
/3 = 2/3, a spin factor of a = 4 and a core radius of 90 
pc. For a more detailed des cription of the hot ga s profile 
and parameters we refer to iMoster et al.l (|2011f ) . These 
sub-halos are expected to have a very low amount of ga s 
(if any), due to re- ionization (e.g. lOkamoto et al.|[2008l ). 


* Actually since we started from a lO"'''J\/0 halo, our halo could 
possibly contain some stars, but in this experiment we decided to 
neglect this component, since it won't change the overall picture. 



Fig. 3. — Effect of the interaction between a DM sub-halo con- 
taining 3.2% of gas by mass(DM mass = 10* Mq). 



Moreover we are interested in "dark" satellites, so we 
do not want the gas in the sub-halos to form stars, 
hence we choose a very low gas fraction (i.e. fraction 
of mass in gas relative to the DM mass) in these sub- 
halos: Mgas/MoM ~ 0.03. We consider only an hot gas 
component; as halos containing cold gas (and possibly 
stars), will be directly detectable; while, in this work, we 
are interested in testing the effects of the more numerous, 
undetected, dark satellite population, which is predicted 
by the cold dark matter model. We then place this new 
sub-halo on the exact same orbit of the pure dark matter 
experiment and let it interact with the gaseous disk as 
described in the previous section. 

The hot gas in the halo is able to displace the gas in 
the disk and produce a low density region (a hole), as- 
sociated with an expanding high density shock wave as 
shown in Figure [3l For our chosen parameters the hole 
has a diameter of about 1.5 kpc and lasts for about 60 
Myr. The lifetime of the holes is decided by two fac- 
tors, the pressure gradient and the differential rotation. 
Our simulations suggest that the pressure gradient causes 
the hole to be filled by gas long before it gets destroyed 
by the differential rotation, although differential rotation 
becomes important for very large holes. For typical hole 
sizes of 1 — 2 kpc we find an average hole lifetime of 
about 50-60 Myr; in agreement with the dynamical ages 
estimated bv lWaher fc Brinkd (|1999[ ). 

Another interesting feature of this simulation is the en- 
hanced star formation on the rim of the hole, as shown 
in Figure 01 New stars have preferentially formed on 
the high density edges of the hole, while the hole itself 
contains very few new stars. This result is in agree- 
ment with the ages of the st ars associated with th e HI 
holes in IC2574 as derived bv iPasquali et all (J2008I ). via 
comparison of observed UBV colors with those predicted 
by synthetic stellar populations (for more details see 
lPasauahet"aI][2008[ ). 

We then run a series of simulations for different per- 
centages of gas fraction of the DM sub-halos, keeping 
all other simulation parameters constant. The effect of 
a pure DM sub-halo and a sub-halo containing gas are 
inherently different. Figure [S] shows the relative surface 
density profiles of the gas along a strip of width 1 kpc 
joining the center of the galaxy and the center of the hole, 
for different sub-halo gas fractions (including the case of 
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Fig. 4. — Star formation rate, showing clearly the ring of new 
stars formed along the rim of the hole. The units of SFR surface 
density in this plot are in 10 A'lQkpc~ . 
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Fig. 5. — The effect of the interaction between DM sub-halos 
with different gas fractions and the disk The curves denotes the 
relative surface density variation with respect to the isolated disk, 
upon the interaction with a halo with 3.2% gas (red), 1.9% gas 
(green), 0.8% gas (blue) and 0% gas (black). The arrow indicates 
to point of impact of the halos. (DM mass = 10* Mq) 



no gas). A pure DM sub-halo (in black) produces a peak 
of increased density while a gaseous halo produces a low 
density region surrounded by high density wave. This 
density wave increases the density in the rim of the hole 
above the star formation density threshold (ptii), thus 
triggering star formation all along the rim of the hole. 
As expected, the size of the hole and the extent to which 
the disk is perturbed depends on the gas fraction of the 
infalling DM sub-halo. 

The impact with a sub-halo will create a net veloc- 
ity gradient along the sub-halo trajectory in the vertical 
component of the H I gas. This signature could in prin- 
ciple disentangle different formation scenarios for the HI 
holes. Figure [6] shows the velocity in the z-component 
(perpendicular to the disk plane - Vz^disk which equals 



Distance from center (kpc) 

Fig. 6. — Variation of V^.disk with initial velocity. The red curve 
traces the change in the velocity profile in the z direction for a halo 
which has 0.8% gas and velocity of v~ = 150 km s~^ and the black 
line for Vz = 300 km s~^, the green curve is for a gas fraction of 
3.2% and velocity of 150 km s~^. 



zero in the unperturbed case) along the same strip as 
mentioned above. Here the infalling halos has two dif- 
ferent gas fractions, 3.2% (green line) & 0.8% (red line) 
with a velocity of Vz ~ 150 km s^^ and a halo with a 
gas fraction of 0.8% with a veloocity of Vz = 300 km 
s~^. A higher sub-halo initial velocity and/or a larger 
gas fraction imparts more energy to the disk in the im- 
pact direction and hence produces a stronger signature in 
Vz.disk- This plot shows that the effect is small and would 
be very hard to detect in real observations as the effect is 
of the same order as the typical velo city resolution of H I 
observations of nearby galaxies (e.g. IWalter et al.ll2008[ ) 

3.3. Cosmological runs 

Our previous simplified setups have shown that dark 
matter sub-halos that contain a small fraction of gas can 
in principle create holes in the HI distribution that re- 
semble the observed ones. The next question to answer 
now is whether there are enough DM satellites to repro- 
duce the number of observed holes in a galaxy similar to 
IC2574. 

To answer this question we turn to cosmological N- 
body simulations which give us the mass, size, veloc- 
ity and orbital parameters of sub-halos of the reference 
galaxy IC2574. We use the publicly available database 
of satellite distribution in the Via Lactea II (VLII) 
simulation (iDiemand et al.l [20071 2008). Unfortunately 
the mass of the parent halo in the VLII simulation 
(2 X 10^'^Mq) is higher than the estimated mass of the 
halo of our reference galaxy IC2574 (IO^^Mq). We thus 
decide to scale down the properties of the sub-halos us- 
ing the following simple dynamically motivated scaling 
relations: 

M' = M/10 (1) 
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(2) 
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Fig. 7. — The velocity as a function of distance for the M12 sub- 
halos(black open squares) after being scaled down and Mil ones 
(red triangles). 
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In order to test that we have not introduced any biases 
in the satellite properties by rescaling them, we run an 
additional N-body simulation of a M ~ IQUMo dark 
matter halo, at a resolution lower than the VLII. We 
have selected the candidate halo from an existing low res- 
olution d ark matter simulation (350^ particles within 90 
Mpc, see iNeistein et all (|2010( )) and re-simulated them 
at higher resolution using the volume rcnormalization 
technique (jKatz fc White! I1993D . The total number of 
particles within the virial radius at z = is « 3.3 x 10^, 
which gives a mass per particle of 3.05 x W^Mq. We will 
refer to this simulation as Mil, while we will use Ml 2 
for the rescaled version of the VLII simulation. Figure 
[7| shows that sub-halos of both runs occupy the same 
phase-space region, which confirms that our simple scal- 
ing is valid as first order approximation Q We use both 
satellite distributions (form Mil and M12) in our cos- 
mologically motivated tests. 

These N-body simulations predict that a large num- 
ber of sub-halos are present around a typical galaxy of 
lOiiA/0. In fact we have 798 sub-halos from the M12 
simulation and 468 from the Mil simulation, with masses 
larger than IO^Mq, our resolution limit in the N-body 
simulations. Not all of these halos pass through the 
disk. Most of them have nearly circular orbits. Only 
a few of them have orbits which take them close enough 
to the center of the halo, and hence interact with the 
disk, these are the halos we arc interested in. We find 
interacting halos by integrating orbits of all sub-halos 
in a static potential for 1 Gyr. This static potential is 
parametrized with an NFW profile, with the same virial 
radius and concentration of the original halo extracted 

^ The shelf of high velocity subhalos at a given distance is due 
to the higher halo statistics in M12, that allows a better sampling 
of the tails (positive and negative) of the velocity distribution. 



TABLE 1 
Number of halos and their properties 

Mass(M0) Nh - M12 Nh - A/11 ^o Ndm N^a 



5.56 X 106 


26 


10 


44.16 


556 


446 


1.01 X lO'' 


7 





41.65 


1010 


799 


1.46 X lO'^ 


3 


2 


40.16 


1463 


1158 


1.91 X 10^ 








- 


- 


- 


2.37 X 10^ 


1 


1 


38.34 


2371 


1876 



Note. — The first column denotes M200 in 
units of solar mass, the second, number of halos 
in M12 run, the third, number of halos in Mil 
run, the fourth, concentration parameter, the 
fifth, number of DM particles used to sample 
the halo and sixth, number of gas particles used 
to sample the halo 

from the Nbody simulation at z = 0. 0. A interaction is 
defined as the passage of the sub- halo within a cylinder of 
radius R = i/xM-y^ < 8 kpc and a height —2<h^<2 
centered at the center of the halo. We obtain a sample of 
37 dark matter sub-halos for M12 and 13 for Mil, where 
the lower number of satellites in the Mil run is due to 
the its lower resolution compared to the VLII simulation. 



In both the M12 and Mil runs DM sub-halo masses 
range between 1.1 x 10*^ A/© - 1.3 x IO'^Mq. Wc divide 
these halos into five halo mass bins. The number of par- 
ticles in each mass bin and their properties are given in 
Table [1] All halos in this run have a gas fraction 5% by 
mass. We put this relatively large amount of gas in these 
halos so as to get an upper limit to the number of holes 
which can be formed. 

For each of the mass bins we generate a synthetic halo, 
using a recipe, which is slightly different from the single 
halo simulation presented in the previous section. We 
create a spherical NFW halo with a density contrast of 
A = 1000 with respect to the critical density of the uni- 
verse, as is typical for sub-halos (the density contrast 
for isolated virialized halos is normally assumed to be 
Amr ~ 100). Each of these halos contains a gaseous halo 
in hydrostatic equilibrium with the DM potential. Fi- 
nally we place back our synthetic halos around the pri- 
mary galaxy at the same position and with the same 
velocity as obtained from the N-body simulations. 

A combination of high mass and high velocity is needed 
to produce holes as shown in the previous section. Due to 
the different orbital parameters such as angle of impact, 
velocity and mass of the halo, the lifetime of the holes is 
more varied, ranging from a mere 10 Myr to as much as 
70 Myr. Figure [S] shows the surface density map of the 
gaseous disk after 0.53 Gyr, a circular hole of about 1 kpc 
can be seen in lower left corner of the simulation. The 
structure seen at the center of the galaxy is the result of 
local instabilities and star formation. 

^ We use as starting point the Nbody results at Zi = 0. Our 
integration time of 1 Gyr would in principle require Zi = 0.079. 
Such a tiny difference should not affect the sub-halo mass function 
on the small scales as considered in this work (M Ri lO''). 
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Fig. 8. — Surface density map of the gaseous disk at a particular 
snapshot of the cosmological run with no winds where the hole 
is clearly seen at (-4,-4). The units of density are the same as 
mentioned in previous surface density plots. 



We now quantify how many holes are produced with 
respect to observations. We focus on large holes (i? > 1 
kpc) , since smaller features can be more easily explained 
by supernovae explosions, stellar winds or a combination 
of both. 

In both our simulation setups (Mil and M12) we ob- 
serve only 3 — 4 events in the 1 Gyr running time (Figure 
[9|) . with no more than one hole in a particular snapshot 
(Figure HI). On the other hand IWalter fc Brinkg (|1999f ) 
observed a total of 8 holes with radii greater that 1 kpc 
in a single "snapshot" of IC2574 (dotted line in Figure 
[9]). The number of holes formed in the Mil run is lower 
primarily due to lack of resolution (when compared to 
M12). This results in an under estimation of the total 
number of low mass satellites. 

It emerges clearly that while DM sub-halos are in prin- 
ciple able to create holes in the HI gas, when a cosmolog- 
ically motivated sub-halo distribution is used the number 
of predicted halos do not match the observed one. This 
leads us to conclude that the sub-halo - disk interaction 
is not the principal mechanism to create the observed HI 
holes. 

We note that the number of holes shown in Figure ^ 
is just an upper limit. So far we have assumed that DM 
sub-halos will keep their gas (needed to displace the gas 
in the H I disc) all the way down to the galaxy center, 
completely neglecting the effect of ram pressure. The 
ram pressure against the hot gas surrounding the galaxy 
could partially or totally remove the baryonic content of 
the sub-halo and make it even more inefficient in produc- 
ing holes. In order to test this scenario we run an addi- 
tional simulation, in which we turn on isotropic galactic 
winds while running the central galaxy in isolation for 1 
Gyr to stabilize it (see section [2]). These winds preferen- 
tially expel gas perpendicular to the plane of the disk as 
it faces the least hindrance in this direction. This cre- 
ates shells of gas surrounding the galactic plane. This 
gas (slightly colder than the gas inside the sub-halos) is 
efficient is stripping all the gas from the sub-halos. 

When this new initial condition (central galaxy run 
with winds on) is evolved together with sub-halos 1 Gyr, 
it produces virtually no holes. This experiment shows 
that it is indeed very difficult for the sub-halos to retain 
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Fig. 9. — Plot showing the cumulative number of holes formed in 
a span of 1 Gyr in the Mil and M12 simulations. The horizontal 
black dotte d line is the number of h olesfd > 1 kpc ) in one snapshot 
of IC 2574 llWalter fc Brinks|[T999h . The lower number of holes in 
the Mil run is, primarily due to its low resolution (compared to 
M12). 



a sufficient amount of gas to perturb the central HI disc. 
This points to a different process being responsible of the 
production of large holes in the HI distribution. 

4. CONCLUSIONS & DISCUSSION 

Atomic hydrogen (HI) observations of nearby galaxies 
reveal complex gas distributions. In many systems, the 
neutral interstellar medium (ISM) contains holes, shells, 
and/or cavities. In order to understand the origin of 
these features, we numerically investigate how the impact 
of dark matter sub-halos orbiting a gas-rich disk galaxy 
embedded in a massive dark matter halo influences the 
dynamical evolution of the HI gas disk of the galaxy. We 
mainly focus our attention on the creation of large holes 
(i? > 1 kpc) in the HI distribution, commonly found in 
observations of nearby galaxies. 

We create a central galaxy resembling the properties 
of the well-studied dwarf galaxy IC2475 and we bombard 
its gaseous disk with dark matter dominated satellites. 
The dark matter, stellar and gas components of both the 
primary galaxy and the satellites are all live (i.e. made 
of particles) and a full hydro-dynamical code (gadget2) 
including cooling, star formation and feedback is used for 
these merger simulations. 

Our experiment shows that a pure dark matter sub- 
halo (M = IO^Mq) is not able to displace the gas in the 
disk, instead due to gravitational focusing it gives rise to 
a localized high density region. This kind of dark impact 
might be able to induce star formation in the outer part 
of the extended gas disk, but cannot create large holes. 
This result was not unexpected: DM particles interact 
only through gravity and do not have the necessary con- 
tact force to push the gas away from the disk. 

These high density tidal imprints can be characterized 
and studied in order to obtain the properties of the im- 
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actin g dark sub- halo as done by I Chang fc Chakrabartil 
201 If ) especially for relatively massive dark satellites 
(1:100). 

To produce holes we need particles which have contact 
forces, hence we assume a small amount of gas to be 
present in the dark matter sub-halo. Even a gas fraction 
as low as 0.8% in a halo of mass IO^Mq with a velocity 
of 150 km/s is able to produce a detectable low density 
region. On average holes have a lifetime of about 20 — 60 
Myr, depending on the halo density, gas mass and impact 
velocity. 

We then use satellite properties directly extracted 
from high resolution Nbody cosmological simulation (Via 
Lactea II, Diemand et al. 2008), to check how many holes 
are predicted in the commonly assumed Cold Dark Mat- 
ter model. 

These cosmological motivated runs show that sub- 
halos with a relatively high gas fraction (5%), arc able 
to produce a total of about 3-4 large holes (i? > 1 kpc) 
in an integration time of 1 Gyr. This number is signif- 
icantly lower than the number of observed holes of the 
same size in IC2475 galaxy. If the effect of ram pressure 
is taken into account, the dark matter satellites tend to 



loose a significant fraction of their gas content, making 
them even less efficient at perturbing the HI disc. 

We conclude that, although DM matter satellites with 
a modest gas content arc in principle able to create holes 
with a radius of order 1 kpc, disk - sub-halo interaction 
are not the primary channel through which these holes 
form in real galaxies. We consider it likely that other 
astrophysical processes like supernova explosions, stel- 
lar winds, high velocity clouds, or a combination of the 
above factors as the main causes for the observed com- 
plex geometry of extended HI disks. 
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